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Special localized wavemodes show up in several physical scenarios including BEC in optical lat- 
tices, nonlinear photonic crystals and systems with strong electron-phonon interaction. These result 
from an underlying nonlinear contribution to the wave equation that is usually assumed to be in- 
stantaneous. Here we demonstrate that the relaxation process of the nonlinearity has a profound 
impact in the wavepacket dynamics and in the formation of localized modes. We illustrate this phe- 
nomenology by considering the one-electron wavepacket spreading in a C60 buckball structure whose 
dynamics is governed by a discrete nonlinear Schroedinger equation with a Debye releixation of the 
nonlinearity. We report the full phase-diagram related to the spacial extension of the asymptotic 
wavepacket and unveil a complex wavepacket dynamical behavior. 
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The discrete nonlinear Schroedinger equation (DNSE) 
has been frequently used to describe a wealthy of nonlin- 
ear physical phenomena associated with wave dynamics. 
The DNSE can display a large class of topologically stable 
solutions such as solitons, vortex rings and breather^"—. 
These solutions have been largely explored in the context 
of Bose-Einstein condensates in optical lattices and light 
propagation in photonic crystals^— . The DNSE also 
captures some features related to the electron-phonon in- 
teraction that can be effectively taken into account by a 
nonlinear term in the electronic Hamiltonian^^. Nonlin- 
earity was shown to breakdown the Anderson localization 
in low-dimensional disordered systems by promoting a 
diffusive- like spread of the electronic wavepacket ^^. 

The self-trapping of the wavefunction is one of the 
most remarkable effects associated with the presence of 
a nonlinear contribution to the discrete Schroedinger 
equation. An initially localized wavepacket spreads in 
the regime of weak nonlinear it ies, but becomes trapped 
around its initial position when the nonlinear coupling is 
above a threshold value^^^^. The self-trapping transi- 
tion has been characterized in one-dimensional and two- 
dimensional systems. The critical nonlinear coupling is 
of the order of the energy bandwidth of the correspond- 
ing linear system for an initial wavepacket fully local- 
ized in a single site. In two-dimensional lattices it in- 
creases as the initial wavepacket becomes wider ^^. Re- 
cently, the influence of the nonlinear response time on 
the self-trapping transition has been investigated in one- 
dimensional systems. By considering a delayed nonlin- 
earity, it has been shown that the critical nonlinearity de- 
pends non-monotonically on the response time^^. When 
the relaxation of the nonlinearity is described by a De- 
bye process, the typical response time drastically affects 
the wavepacket dynamics. Much weaker nonlinear it ies 
are needed to trap the wavepacket ^^. It focalizes after an 
initial spread exhibiting pronounced finite-size effects in- 
cluding wavepacket fragmentation. The relaxation of the 
nonlinearity also leads to the wavepacket relocalization 
in disordered systems— i^. 



The wavepacket dynamics in nonlinear honeycomb and 
some closely related lattices has been a subject of in- 
tense recent investigations. These have been mainly 
motivated by the ability to build up optical lattices 
and photonic crystals with these topologies, which al- 
lows the experimental observation of several theoret- 
ically predicted phenomena, such as topology-induced 
bistability^, discrete localized modes and breather o^^i^^ , 
conical diffractio n^^i^^ , self- trapping of vortices, gap 
vortices and gap solitons^S.. Nonlinearity is also ex- 
pected to play a significant role in the dynamics of elec- 
tronic wavepackets in carbon-based structures^i^. Re- 
cent experiments have probed the electron-phonon in- 
teraction in graphene which can accounts for an effec- 
tive nonlinearity^^. Low-temperature scanning tunneling 
spectroscopy experiments have been recently used to map 
the wavefunction in graphene quantum dots^^, which 
open the possibility to experimentally observe many the- 
oretically predicted nonlinear effects in low-dimensional 
carbon-based structures. However, most of the previous 
theoretical studies of nonlinear wave dynamics in two- 
dimensional lattices do not include finite-size effects and 
the relaxation of the nonlinearity which might play im- 
portant roles in nano-scaled structures. 

In the present work, we show that both finite-size and 
finite response-time effects have a strong impact on the 
wavepacket spreading dynamics. We will follow the time- 
evolution of an initially fully localized wavepacket whose 
propagation is governed by a DNSE with Debye relax- 
ation of the nonlinearity. To illustrate the phenomenol- 
ogy, we consider the wavepacket evolution over a C60 
buckball. In this topology, all sites are equivalent and 
the absence of borders avoids the superposition of surface 
effects to the finite-size ones. We will report the phase 
diagram in the two dimensional parameter space (nonlin- 
ear couphng X and response time r) which unfolds a quite 
complex structure with three distinct localized phases. 

A nonlinear contribution for the one-electron 
Schroedinger equation results from the underlying 
electron-phonon interaction——. Within the adiabatic 



approximation, such nonlinear term can be considered as 
instantaneous, with the on-site potential at a given time 
depending on the local electronic density at this proper 
time. Going beyond the adiabatic approximation, the 
effective nonlinear term has been shown to obey a Debye 
relaxation process^^. Within a tight-binding approach 
and considering a localized orbitals basis, the dynamics 
of a one-electron wavepacket can be described by a 
DNSE with a relaxing nonlinearity, written as 

m 

Xn{t) = --[X„(t) + x|*„(t)P] , (2) 

where we used h = 1. ^n(^) is coefficient of the wave 
vector expanded in the localized orbitals basis |^(t)) = 
5]]^ ^n(^)l^)- The relaxing nonlinearity X{t) is consid- 
ered to have a typical response time r and strength co- 
efficient X- The sum is taken over first-neighbor pairs 
of sites (n, m) that are coupled by a hopping ampli- 
tude Vnrri' ^n's represent the on-site energies. In the 
limit of r -^ 0, the nonlinearity becomes instantaneous 
Xn{t) = -x|^n(^)P and the usual DNSE is recovered. 

In what follows, we will solve the above nonlinear dy- 
namical equations for an initially localized wavepacket 
^n(^ = 0) = ^n,i- We will consider a C60 buckball topol- 
ogy (see fig. 1) which is composed of 20 hexagons and 12 
pentagons. Each of the 60 equivalent sites have three 
bonds. Two of them are between a hexagon and a pen- 
tagon and the remaining one between two hexagons. In 
the real C60 molecule, the hexagon- hexagon bonds are 
slightly shorter than the hexagon-pentagon bonds. This 
feature ultimately leads to a small difference in the corre- 
sponding hopping amplitudes which we will not consider 
in the present work. We will take Vnm = V for every 
bond and will work in units of V = 1. Further, as all 
sites are equivalent, we can also consider £n = s = 
without any loss of generality. The dynamical equations 
were solved using a standard eighth-order Range-Kutta 
algorithm and the norm conservation was followed with 
a precision of 10~^ to ensure the numerical accuracy. 

In order to characterize the spacial extension of the 
wavepacket, we computed the participation function de- 
fined as 
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which gives a measure of the fraction of the lattice sites 
at which the wavepacket is spread at time t. It becomes 
equal to N for an uniformly distributed wavepacket. For 
a fully localized state P — 1. 

We followed the time-evolution of the wavepacket until 
it reached a statistically stationary regime. We explored 
the two-dimensional parameter space (r, x) in order to 
build up the full phase-diagram regarding the spacial 
extension of the asymptotic wavepacket. Four distinct 
phases were identified, as reported in figure 2. There is a 







FIG. 1: (Color online) Schematic representation of the C60 
buckball. Each site has two hexagon- hexagon bonds and one 
hexagon-pentagon bond. Ci represents the site at which the 
wavepacket is fully localized at t = 0. C2 is the site diamet- 
rically opposite to the initial position. C3 stands for the site 
targeted in figure 3d. 



single extended phase, on which the wavepacket spreads 
uniformly over all sites of the buckball, and three dis- 
tinct localized phases. The extended phase prevails in the 
regime of weak nonlinear it ies. However, narrow branches 
of the extended phase are found between localized phases. 
While these extended branches die out when the relax- 
ation time increases, the main extended phase becomes 
wider. 

The predominant localized phase Li corresponds to an 
asymptotic wavepacket that results localized mostly at 
the initial position (site labelled Ci in fig. 1). In this 
region, we found the participation number to become 
slightly above unit, signalling a quite localized state, with 
the probability density concentrated at the initial posi- 
tion. There are two main regions of the parameter space 
on which the wavepacket becomes trapped on the initial 
position. A second localized phase L2 appears in be- 
tween these two Li regions. In this new localized phase, 
the participation function is of the same order as in the 
localized phase Li. However, the probability density is 
centered at the site that is diametrically opposite to the 
initial position (labeled as C2 in fig. 1). The two branches 
of the extended phase are in between these two localized 
phases. A third localized phase L3 occurs in between Li 
and L2. In this narrow phase, the wavepacket was found 
to become localized around other positions besides the 
initial and the diametrically opposite sites. It is interest- 
ing to call attention to the fact that a complex sequence 
of localization-delocalization transitions takes place at 
small values of the relaxation time when one crosses the 
regime of intermediate nonlinear strengths. 

In order to have a clearer view of the wavepacket 
dynamics on these distinct phases, we plot the time- 
evolution of the participation function and the relevant 
occupation probabilities for some representative values of 
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FIG. 2: (Color online) Phase diagram concerning the spa- 
cial extension of the asymptotic wavepacket in the (r, x) pa- 
rameter state. The extended phase with a uniformly dis- 
tributed wavepacket occurs at small nonlinearities and two 
branches. Li corresponds to a phase with the wavepacket 
localized around its initial position. L2 is the phase where 
the localization is concenterated in the site diametrically op- 
posite to the initial position. In the narrow L3 phase, the 
wavepacket localizes in other sites. A complex sequence of 
localization-delocalization transitions is observed at small r 
values and intermediate nonlinear strengths. 



the physical parameters, as shown in fig. 3. The dynam- 
ics in the extended phase is illustrated in fig 3a for three 
set of parameter values which cover the two branches 
and the main extended region. The normalized partici- 
pation function converges to P/N = 1 indicating that the 
asymptotic wavepacket is uniformly distributed over all 
buckball sites. The convergence to this state is achieved 
after an oscillatory transient regime, thus indicating that 
the expanding wavepacket has a breathing character. 

In fig. 3b we depict the typical wavepacket dynamics 
in the Li localized phase. We show the time-evolution of 
both the participation number and the occupancy prob- 
ability of the initial site |^i(t)p. In the initial tran- 
sient, the wavepacket spreads in an oscillatory fashion 
and occupies a finite fraction of the buckball sites. After 
this transient, the wavepacket self- focuses and relocalizes. 
The increase of |^i(t)p indicates that the wavepacket re- 
localizes around its initial position. This relocalized state 
was found to be stable over a very long run. A very sim- 
ilar phenomenon occurs in the localized phases L2 and 
I/3, as illustrated in fig. 3c and fig. 3d respectively. How- 
ever, the self- focusing drives the wavepacket towards the 
site diametrically opposite to its initial position in phase 
1/2. In the localized L3 phase, the wavepacket spreads 
almost uniformly over the lattice before starting to self- 
focuses. In this case, the position of the asymptotic lo- 
calized state is very sensitive to values of the physical pa- 
rameters. For the values used in fig. 3d, the wavepacket 
relocalizes around the site labeled as C3 in fig. 1. 

Before finishing, we closely analyze the wavepacket 
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FIG. 3: (Colour online) Time-evolution of the normalized par- 
ticpation number P/N and the relevant occupancy probabili- 
ties |^i(t)|^ for representative values of the physical parame- 
ters, (a) Extended phase: x — 5-5 and r = 0.22, 0.45 and 1.50 
(from right to left) being respectively on the first, second and 
third extended regions. Notice that the expanding wavepacket 
develops breathing oscillations, (b) Localized Li phase: after 
the initial expansion, the wavepacket self- focuses and relocal- 
izes around its initial position, (c) Localized L2 phase: in this 
case, the wavepacket relocalization occurs around the site dia- 
metrically opposite to its initial position. Localized L3 phase: 
here the relocalization occurs around the site labbelled as C3 
in fig. 1. 
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FIG. 4: (Color online) (a) Asymptotic normalized participa- 
tion function P/N an occupancy probabilities |^i(t)|^ at the 
initial (i = 1) and diametrically opposite {i = 2) sites as a 
function of the nonlinear coupling for r = 0.0625. The system 
depicts a complex sequence of localization and delocalization 
transitions, (b) The time-evolution of the normalized par- 
ticipation function for r = 0.0625 and two slightly different 
nonlinear strengths. Notice a dynamical transition that leads 
to distinct phases. 



dynamics in the parameter region where a complex se- 
quence of localization-delocalization transitions occurs 



within the main extended region (see fig. 4). We firstly 
plot the asymptotic normalized participation function 
and occupancy probabilities of the initial and diamet- 
rically opposite sites as a function of the nonlinear cou- 
pling and a fixed relaxation time r = 0.0625, as shown 
in fig. 4a. P/N = 1 in the extended phase, as ex- 
pected for a uniformly delocalized wavepacket. However, 
before the ultimate localization transition that traps 
the wavepacket around its initial position, a complex 
sequence of localization-delocalization transitions takes 
place. The position of the localized wavepacket appears 
to randomly jump between the initial and diametrically 
opposite sites. For other small values of the relaxation 
time, we verified that the wavepacket can eventually re- 
localizes around other sites. The high sensitivity of the 
wavepacket dynamics on the precise values of the param- 
eters in this region is illustrated in fig. 4b. One observes 
that, while the orbits of two cases with slightly differ- 
ent parameters remain quite close to each other during 
some time, a dynamical transition takes them apart. It 
is worth to recall that a similar dynamical transition has 
been previously reported to occur in a two-level system 
driven by a relaxation process of the nonlinearity^^. 

In conclusion, we showed that the wavepacket spread- 
ing dynamics displays a rich phenomenology including 
self-trapping, defocusing and focusing when the relax- 
ation process of an effective third-order nonlinearity is 
taken into account. A C60 buckball structure was used 
to illustrate this phenomenology. By considering a fully 
localized initial state, we provided the full phase diagram 
for the spacial extension of the long-time wavepacket. 
The phase diagram was shown to have a complex struc- 
ture unveiling the presence of three distinct phases of lo- 
calized wavepacket s, besides the extended phase in which 
the wavepacket spreads uniformly over all buckball sites. 
The extended phase predominates at small nonlineari- 
ties. In the localized phases, the wavepacket was shown 
to relocalize after an initial spreading over a finite frac- 



tion of the network sites. Besides the usual relocalization 
of the wavepacket around its initial position, focalization 
on the diametrically opposite site takes place in a large 
stripe of the phase diagram. A narrow phase where the 
wavepacket relocalizes around other sites was also iden- 
tified. A complex sequence of localization-delocalization 
transitions was shown to occur at small relaxation times 
as a function of the nonlinear coupling. These were as- 
sociated to a dynamical transition which is strongly sen- 
sitive to the parameter values in this region. A stability 
analysis of the stationary states would be in order to 
shed light in the underlying mechanism that drives these 
transitions. 

The main aspects of the here reported phenomenology 
are expected to show up in general nonhnear physical sys- 
tems where the wavepacket dynamics is driven by a re- 
laxation process of an underlying nonlinearity. These in- 
clude electronic states in nano-sized clusters with strong 
electron-phonon coupling, BEG in optical lattices and 
light propagation in nonlinear photonic crystals. The 
present work leaves open the quest of investigating the 
possible influence of open boundaries and saturation ef- 
fects that might also be relevant in optical and BEG ex- 
periments in finite planar lattices. With the recent ad- 
vances of experimental techniques that allows to map a 
diversity of localized wave modes, further theoretical ef- 
forts would be in order to evaluate the relative relevance 
of these ingredients to the wavepacket dynamics in the 
presence of a relaxing nonlinearity. 
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